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We analyze a trapping reaction with a single penetrable trap, in a one dimensional lattice, where 
both species (particles and trap) are mobile and have a drift velocity. We obtain the density as seen 
from a reference system attached to the trap and from the laboratory frame. In addition we study 
the nearest neighbor distance to the trap. We exploit a stochastic model previously developed, and 
compare the results with numerical simulations, resulting in an excellent agreement. 
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> ■ I. INTRODUCTION 

o : 

I With the aim of obtaining a description of difFusion-hmited reactions several generalizations of the classical Smolu- 
. chowski model have been proposed A stochastic model (SM) coming from the theory of nuclear reactors 
^-H '■ has been adapted to model these reactions [Hjs). This framework results to be very adequate to describe different 
situations, making it possible to obtain results in excellent agreement with simulations j^-p^, particularly when the 
model equation can be exactly solved. 

The trapping reaction A + B B with a single trap and a uniform initial density of particles in a one-dimensional 
lattice is one of the simplest cases describing a certain kind of chemical reaction. In the present work we analyze this 
situation within the framework of the SM, when both species are mobile and have a drift. 
T-H ■ Diffusion-controlled reactions in presence of a biasing field are of great interest in many areas. Let us mention 
'""I ' a few instances: (a) Due to the gravity field, all diffusion processes that take place over the earth are affected in 
some quantity. For exampl e, g ravity is responsible for the drift of water molecules diffusing underground Jl5| , p^ . 
(b)In positron tomography [l^ the addition of a strong field perpendicular to the surface leads to greater sample 
penetration, (c) In problems of electromigration jl^ the model discussed here could be used to analyze the effect of 
a trapping centers on carrier concentration. 
^ ' A situation related to the presented here was studied in Ref. where the case of a stationary and impenetrable 
trap was analyzed. In the mentioned work a diffusion equation with adequate boundary conditions is the base of 
(~| calculations. We want to remark that in this paper we model a standard trap, where the particles can eventually pass 
O from left to right of it (or vice versa) without being trapped. As we will see, while calculations where an impenetrable 
^ trap is involved are naturally treated with a diffusion equation, standard trap problems are naturally treated within 
^ , the SM. This difference is not relevant for particles without a bias in its diffusive motion, at least for the density 
calculations. 

The organization of the paper is as follows. In the next section we briefly describe the model and discuss the 
' differences between penetrable and impenetrable traps. In Sec. Ill we calculate the density as seen from a reference 
■ ■ ' system attached to the (mobile) trap and discuss some symmetries. In the next section we make a similar analysis 
for the standard density, i. e. from the laboratory reference system, and calculate the number of particles absorbed 
by the trap. In Sec. V we present the results for the nearest neighbor distance. Finally we draw our conclusions. 



II. THE MODEL 



The model considers two species of particles, A and B, both mobile and independent, and having a given reaction 
probability when they meet. Here we will study a one-dimensional trapping reaction in a system of A particles 
uniformly distributed and a single trap i?, with both species having a biased diffusive transport. 
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The model equation for the evolution for N{x,t), the density of the A particles for a given realization of the trap 
trajectory e{t), is the following 

^7V(x, t) = Da^N[x, t) - VA^Nix, t) - 75 (x - e(i)) N[x, t), (1) 

where 7 is a constant measuring the reaction probability (7 — > 00 for perfect reactions), Da,b is the diffusivity and 
va,b characterizes the bias of the A ov B particles. In general we are interested in n{x^t) = {N(x,t)), the density 
averaged over realizations of e{t) in the laboratory frame, or the similar (averaged) density n{x, t) seen from a reference 
system attached to the trap. 

As was pointed in Ref. |^ the SM equations given in Ref. [Q are valid using the propagators corresponding to the 
particle transport process. In our case we have for A and B particles 
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respectively. 

We will consider for a moment that the trap is fixed at the origin and the A particles are not biased. Then Eq. (|^) 
(in this case N{x,t) = n{x,t) = n{x,t)) becomes 

d 

— n(a;, t) = DA-^n{x, t) - ^5 (x) n{x, t). (4) 

The resulting equation is equivalent to a diffusion equation except at the origin. It is not difficult to obtain the jump 
of the derivative in the density which due to the symmetry results in the following boundary condition (for x > Q) 

d 7 

—n{x,t)\^^o = ^jj-n{0,t). (5) 

Hence the diffusion equation with this boundary condition is equivalent to Eq. (^). The same boundary condition could 
be obtained requiring that the diffusive current at the trap position (left minus right contributions) is proportional 
to the density, which is the usual way to obtain the so called radiative boundary condition for an imperfect trap. 
However, this is true due to the symmetry, which requires that the current at both sides of the trap are of equal 
magnitude. When coming to our bias problem (but with the trap immobile) the particles have a preferred direction of 
motion and the symmetry argument is not valid. On one hand, Eq. (0) models a penetrable trap but, although we can 
write an equation for the jump of the derivative we can not to write a equation for each lateral derivative separately. 
On the other hand, a current proportional to the density (from one particular direction) models an impenetrable 
trap. In this way we can see that each model is the natural way to treat each different situation. It is intuitive that 
for perfect trapping (7 — > 00) both models give the same density. A final remark is that the coefficient relating the 
diffusive current and the density at the trap position was taken as 7. This election make both models compatible. 
We must change 7 by j/{2Da) in the equations of Ref. pof in order to make a direct comparison with results given 
in this work. 

A detailed discussion for unbiased particles (and fixed trap) of penetrable and impenetrable traps was done by 
Taitelbaum Note that the analysis done here about impenetrable traps is only for a fixed trap, the case considered 
in Ref. . The treatment for a mobile impenetrable trap is a non trivial generalization, and is out of the scope of 
this work. 



III. DENSITY IN THE TRAP FRAME 

In complete analogy with the case without bias we change the reference system to one fixed to the trap. 
After doing the average it results that the density is similar to the case with a fixed trap, but changing Da by 
Di = Da + Db and va hy Vr = va — vb- The resulting equation is solved by Laplace transform techniques, yielding 
after some calculations the expression 

n{x,t) l^vAx-\x\)/(2D,) i^vAx+\x\)K2D,) f\x\+Vrt\ 

— 1 — erfc — — eric — (6) 
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valid for any value of the parameters except if 
obtained performing a specular reflection) 



7. For Vr = J the expression is (expression for Vr = —7 can be 



"(^'^) ^ 1 _ levAx-\x\)/(2D,)^^^^ nx\-Vrt\ 
no 4 V V^D^i J 



(7) 



We can see that n depends on the drift velocities only through the relative velocity u^, and when both species have the 
same drift we recover the unbiased case. [ p"3|j2l| ] In addition, a change in the sign of Vr produce a specular reflection 
in the density proflle. Hence, without loss of generality in the analysis we often set Vr > 0. Note that in Eq. (^) the 
variables and parameters come only in three combinations: x/y/Dit, Vr\/t/Di and 7/fr- 

Figure 1 shows a typical n at different times and four sets of simulations corresponding to the same theoretical 
profile as all have the same Di and Vr. The agreement between the SM result and simulations is excellent. 

We can observe that for t 00 the density in the trap position approaches 



no 



7+ l^rl' 



(8) 



which is a positive value (for imperfect trapping), in contrast with the unbiased case {vr = 0) where the particle 
density asymptotically reaches zero at the trap position. [ p"3pl| ] Besides, for Vr > the density reaches a stationary 
profile {n{x,t)/no ~ 1 — exp{vrX / Di) / {1 + fr/7)) in the left semi-axis, i. e. where the relative drift is towards the 
trap, while any point to the right of the trap reaches asymptotically the value given by Eq. (^). These aspects can 
be seen in Fig. 1. 



IV. DENSITY IN THE LABORATORY FRAME 



The Laplace-Fourier transform of the density in the laboratory frame H] is given by 
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Here we can see that Eq. (^ does not change if we interchange Da by Db and va by vb simultaneously. This fact 
shows that the trap and the particles play exactly the same role in the density. However, this characteristic is not 
valid for any magnitude, as is known from the first neighbor distance in the unbiased problem p^ . The equivalence 
between both particle species just mentioned lead to two interesting consequences that we describe in the following. 

The first consequence is a symmetric density if we set Da = Db = D and va ~ —vb = v. This can be verified 
considering that from Eq. (^ results n{k, s) = n{—k, s). To obtain the general expression in the x and t variables is 
very difficult and we present the simple result for perfect absorption 
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This expression has only two variables (any parameter) that are x/y/ADt and v^Jt/(AD). We note that this expression 
approaches asymptotically zero at the origin, in contrast with the unbiased case where the density reaches a finite 
value. |]l3| , pT| In Fig. 2 we can see an example of symmetric profiles at different times. The agreement between theory 
and simulation is good, although there is a systematic difference. However, note that the symmetry is still present in 
simulations. 

The second consequence is that setting Db ~ vb ~ (flxed trap) the density expression is the same as setting 
Da = Va — (fixed particles) exchanging the particle subindexes. Moreover, since for the fixed trap we have the 
trivial equality h{x,t) = n(x,t), and n(x,t) for Db = = is equal to n{—x,t) for Da = va = ^, we have 
h{—x,t) = n{x,t) for fixed particles, which is not evident. A similar property was found in the unbiased case [O. 
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In Fig. 3 we show simulations of fi{x, t) and n(a;, t) for fixed particles, n(x, t) for the fixed trap and the SM result 
corresponding to all of them. We can see an excellent agreement. 

By setting fc = in the second term of Eq. the total number of absorbed particles in the Laplace domain can 
be obtained. Performing the inverse transform we obtain 
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The limit cases 7 cx) and 7 = not included for the sake of brevity, can be extracted from Eq. ( |ll|) taking 
the corresponding limit. In Fig. 4 we show simulations of Nabs for the same sets of parameters than in Fig. 1. 
Simulation data coming from all sets are rather coincident (almost indistinguishable in the figure) and in excellent 
agreement with the SM prediction. If — vb Eq. ( pi] ) coincides with the case without bias, where asymptotically 
Nabs oc t^'^ [pi. However for Vr — VA — vb ^ the asymptotic behavior results 
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This change in the behavior can be understood from Eq. (g), since the density at the trap position reaches a positive 
value. Hence, the absorption rate 771(0, i) becomes constant. The same arguments are applicable to an impenetrable 
trap that has the same qualitative dependence (although the factor is different because n{0,t) is different). 

It is possible to factorize Eq. (|^) (after performing the inverse Laplace transform) in exp{—ikvBt) and an expression 
where va and vb arise only in the Vr combination. This implies that in the x domain the density has the form 
n{x,t) — n*(x — VBt), where n* depends on Vr but neither on va nor vb separately. Due to the identical role played 
in the density expression by both particle species, we can do a similar factorization involving exjp(—ikvAt). 

In Eq. (^ we take the limit 7 — > 00 in order to obtain the mathematically simple expression for perfect reaction 
resulting, after performing the inverse Laplace transform, in 
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Numerically performing the inverse Fourier transform we obtain the density in the x domain. In Fig. 5 we show 
numerical integrations at different times together with the respective simulations. The agreement between both results 
is apparent. 



V. NEAREST NEIGHBOR-DISTANCE 

As was stated in Ref. [|l^, we can use the expression for the PDF calculated in a similar way as for fixed particles 
as an approximation to the problem of both mobile species. This approximation is better for a small Db/Da ratio 
and for short times. We note that we can consider biased particles as having a deterministic drift velocity plus an 
unbiased diffusion. Hence, following Ref. ||l^ the PDF only depends on Vr (at variance with diffusivities where the 
dependence is not trivial). To calculate the PDF of the (one sided) nearest neighbor distance we first compute the 
intermediate quantity Q{x,t) = exp (^— n{u,t)du) , and then the PDF is p{x,t) = —dxQ{x,t) = h{x,t)Q{x,t). 
|p2|-p4| The integral involved in Q can be done analytically, yielding a rather complicated expression. 

In Fig. 6 we show a plot of p corresponding to the right neighbor and compare it with simulations. We can see a good 
agreement. The most noticeable difference arises near the origin. However, the mean distance (x) = p{x,t)xdx — 
/o°° Q{x,t)dx shows the expected systematic differences due to the approximation. This is apparent in Fig. 7 where 
we show (x) coming from a numerical integration. Note that squares and diamonds are almost coincident in Fig. 7 
(because they have the same parameters except the drift velocities but with the same Vr). In addition we make a 
similar calculation and simulations for the left neighbor. In this case the mean distance reaches a stationary value 
that can be calculated within the SM and is 
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where 7(0,2) is the incomplete gamma function |25| (we have the over bar to avoid confusions with the absorption 
constant). We can see in the figure that this quantity does not have a good agreement with simulations. This is 
due to the fact that the model is a continuous one and can not describe well the situation very near the trap, where 
the discreteness of the problem is relevant. This phenomenon is the same as in Fig. 6 (small x), but here it is more 
noticeable because almost all particles in the simulation arc in the same place of the trap or in the nearest neighbor 
site. 

Simulations were performed on a lattice of L sites, with periodic boundary conditions. We have used the same 
algorithm described in jl^], with a small change to allow for the drift of the particles. The connections between the 
simulation and SM parameters are D = ojjAx'^ /2, v ~ {2q — l)ujjAx and 7 = Wr-Acc; where Ax is the jump length, 
ujj^r are the jump and reaction frequencies respectively, and q is the probability of a given jump to the right {q = 0.5 
corresponds to an unbiased particle). All show simulations are the average of 10 000 realizations over a 200 site lattice 
and with an initial density uq — 1. 



VI. CONCLUSIONS 



We have studied a trapping reaction A + B B with a single penetrable trap i? in a one dimensional lattice, where 
both species are mobile and perform a biased diffusion. 

We have seen that the natural tool to perform analytical calculations is the stochastic model (SM) in contrast with 
a diffusion equation that is useful for calculations involving an impenetrable trap. We have calculated the expressions 
of fi{x,t), the distribution of particles in the trap frame, and n(x,t), the usual density in the laboratory frame. 

Analogously to the unbiased case, h{x, t) can be obtained from the fixed trap case doing the exchange Da Di = 
Da + Db and va Vr — va ~ vb in the corresponding expression. In the general expression we note that n(0,t) 
for imperfect reactions reaches a non zero value at long time. This results in a linear asymptotic number of absorbed 
particles, in contrast with the t^/^ behavior of the unbiased problem. 

The expression for n{x, t) has the following symmetry: it is invariant by exchange of diffusivities and velocities 
simultaneously, i. e. Da Db-, Db Da, va vb, vb va- This is valid both for perfect and imperfect 
reactions. Due to this characteristic there is a combination of parameters, namely Da = Db and va — ~vb, which 
results in a symmetric density. Another consequence is that n{x,t) = n(—x,t) for fixed A particles, a property found 
first in the unbiased case. We note that in the symmetric case for perfect reaction the particles are asymptotically 
extinguished, in contrast with the unbiased case where the density reaches a positive value. An explanation of this 
fact is that the trap in average goes to the right and particles located to its right can not pass through it to reach the 
origin (because it is a perfect trap). In addition particles situated to the left have in average a velocity to the left. 
Hence the origin results asymptotically free of particles yielding the observed null density. 

A comparison of the density profile with that obtained for an impenetrable trap, when the trap is fixed (where 
we have available results from Ref. |]l^), shows that both expressions are similar if the absorption rate is high. In 
fact, for perfect absorption both densities become identical. On the contrary, for low absorption an impenetrable trap 
implies that the untrapped particles are accumulated (to the left if va > 0) while they follow their way to the right 
for a penetrable trap. In spite of this, both problems reach a stationary profile to the left of the trap (for va > 0). 

We have studied the nearest neighbor distance in the case of Db =0, and used an approximation valid for small 
Db- We can see that the mean distance, after a short transient, approaches a constant to the left of the trap, while 
it grows to the right of the trap. 

Simulations arc in excellent agreement with the model. As was mentioned in the introduction, this is an expected 
result since the model has been solved exactly. Finally, we want to remark that while in unbiased particle problems 
the SM is a convenient alternative to other technics, in problems of biased particles it becomes the natural way to 
obtain analytical results. 
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FIG. 1. Density of A particles in the reference frame of the trap at different times. Sirrmlations are performed with different 
parameter sets indicated in the figure. All of them have D\ = 0.2 and Vr = 0.16. The reaction is imperfect with 7 = 1. 

FIG. 2. Density of A particles in the symmetric case, as seen from the laboratory coordinate system, for a perfect reaction. 
The simulation parameters are: Da = Db = 0.1 and va = —vb = 0.08. 

FIG. 3. Simulation for the density of fixed A particles in the reference frame of the trap, from the laboratory coordinate 
system, and for the fixed trap. The parameters are: Da or Db = 0.1, \va\ or \vb \ = 0.08 and 7 = 1. 



FIG. 4. Number of absorbed particles. The parameters are the same as those of the Fig. 1 



FIG. 5. Density of A particles, as seen from the laboratory coordinate system, for a perfect reaction. The simulation 
parameters are: Da = 0.2, Db = 0.1, va = 0.16 and vb = —0.08. 



FIG. 6. PDF of the nearest neighbor distance to the trap at t = 200. The parameters, indicated in the figure, are the same 
as in Fig. 1. 



FIG. 7. Mean distance to the nearest neighbor to the trap. The parameters, indicated in the figure, are the same as in Fig. 

1. 
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